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(^ ' Abstract. We have observed the cluster Abell 2218 (z = 0.175) with ISOCAM on board the Infrared Space Observatory using 

^^ . two filters, LW2 and LW3, with reference wavelengths of 6.7 and 14.3 yum, respectively. We detected 76 sources down to 54 

^^ ' and 121 yujy (50% completeness levels) at 6.7 and 14.3 yum, respectively. All these sources have visible optical counterparts. We 

(-H ^ have gathered optical and near-infrared magnitudes for 60 of the 67 non-stellar optical counterparts to the ISOCAM sources, 

O ■ as well as redshifts for 43 of them. We have obtained acceptable and well constrained fits to the observed spectral energy 

distributions (SEDs) of 41 of these sources, using the "GRASIU' models of Silva et al. (1998), and have determined their total 
infrared luminosities {Lju's) and star formation rates (SFRs). 
The SEDs of 20 (out of 27) ISOCAM cluster members are best fit by models with negligible ongoing star formation, and no 

—t I major episode of star formation in the last ~ 1 Gyr. Their SEDs resemble those of 5-10 Gyr old early- type galaxies. A slightly 

higher, but still very mild, star-formation activity is found among the remaining cluster sources, which are mostly spirals. The 
median IR luminosity of the 27 ISOCAM cluster sources is Ljj; = 6 x 10** Lg. The ISOCAM-selected cluster galaxies have 
indistinguishable velocity and spatial distributions from those of the other cluster galaxies, and do not contribute significantly 
to the Butcher-Oemler effect. If A2218 is undergoing a merger, as suggested by some optical and X-ray analyses, then this 

C^ ' merger does not seem to affect the mid-infrared properties of its galaxies. 

The SEDs of most ISOCAM-selected field sources are best fit by models with moderate ongoing star formation, with a sig- 
nificant fraction of their stellar mass formed in the last ~ 1 Gyr Their SEDs resemble those of massive star-forming spirals 
or starburst galaxies, observed close to the maximum of their star formation activity, but not necessarily during the short-lived 
starburst event. The median redshift of these field galaxies is z == 0.6. Their L«'s span almost two orders of magnitudes, 
from ~ 10'" Lq to ~ 10'- Lq, with a median of 1.2 x 10" (eight of the 14 field sources are LIRGs). The SFRs of these 14 
ISOCAM-selected field sources range from 2 to 1 25 Mq yr" ' , with a median value of 22 Mq yr" ' . 
We compare our findings with those obtained in other ISOCAM cluster and field surveys. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual (Abell 2218) - Infrared: galaxies 

1. Introduction 

Clusters of galaxies are not a very common galaxy environ- 
Send offprint requests to: A. Biviano m^nt, but they are crucial for our understanding of galaxy for- 

* Based on observatios with ISO, an ESA project with instruments motion and evolution. In a hierarchical cosmological scenario, 
funded by ESA Member States (especially the PI countries: France, Ae first galaxies to form are those located in high density 
Germany, the Netherlands and the United Kingdom) with the partici- peaks, which then evolve to become galaxy clusters (e.g.. Coles 
pation of ISAS and NASA et al. 1999) . Hence, on average, cluster galaxies are believed to 
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be older than field galaxies (e.g., Kauffmann & Charlot|1998l. 
When a field galaxy is accreted by a cluster, it is likely to speed 
up its evolution rate, as a consequence of the physical processes 
which are switched on in the new, rather extreme, environmen- 
tal conditions, such as ram-pressure, tidal stripping, 'harass- 
ment', and galaxy-galaxy collisions (see, e.g. Abraham et al. 
[19961 Acreman et al.l2003', Gnedin|2^. While the long-term 
effect of these processes is to transform a star-forming galaxy 
into a quiescent one (or to disrupt it entirely, in the case of 
dwarfs), it is yet unclear how this transformation occurs, and, 
in particular, if the affected galaxy undergoes a starburst phase 
during the transformation process (see, e.g., Poggianti 120031 
Okamoto & Nagashima 2003 1. 

A direct observational evidence of the evolution of the 
properties of cluster galaxies with redshift z, is the Butcher- 
Oemler (BO, hereafter) effect , i.e. the excess of blue galaxies 
in distant clusters relative to nearby ones (see, e.g.. Butcher & 
Oemler|T984 Margoniner et al. l20Un Pimbblet"2003 1. The ef- 
fect is partly, yet not entirely, due to field contamination (see, 
however, Andreon et al. 120041 who argue that the BO-efFect 
could be due to observational biases). The BO-efFect is based 
on photometric data. The spectroscopic analogue of the BO- 
efFect is the increase of the fraction of 'k-i-a' galaxies (origi- 
nally called 'E-hA' by Dressier & Gunn 1983 1 with z- The spec- 
tra of k+a galaxies are characterized by strong Balmer lines 
in absorption and the absence of (strong) emission lines (see 
e.g. Dressier et al. T?99 Ellingson et al. 2001 1. These spec- 
troscopic features are usually interpreted as the signatures of 
a post-starburst phase (see, e.g. Poggianti I2003> . Another di- 
rect evidence of the evolution of cluster galaxies is the change 
in the morphological mix of clusters with z (Dressier et al. 
[T9971 Fasano et al. l2000l : clusters at 0.2 < z < 0.6 have a 
significantly lower fraction of SOs (or early-type galaxies, see 
Fabricant et al. ■2000J than nearby clusters. This suggests that 
a transformation of spiral galaxies into SOs took place over the 
last ~ 4 Gyr 

It is yet unclear how the BO-efFect is related to the change 
in the morphological fractions of cluster galaxies. The blue BO 
galaxies seem to be mostly disk galaxies with disturbed mor- 
phologies (Lavery & Henrv .l994j . They could be spirals seen 
before they transform into SOs. If the transformation occurs af- 
ter the gas content of these galaxies is used in one or several 
starbursts, an excess population of post-starburst galaxies is 
naturally expected at an intermediate phase of the transforma- 
tion process. Since starbursts usually occur in very dusty envi- 
ronments, it is important to explore these processes at infrared 
(IR) wavelengths, where the stellar radiation is re-emitted, in 
the presence of dust (see, e.g. Poggianti 2003 1. In particular, it 
was shown by Poggianti et al. C 1999. 1 that a particular spectral 
class of cluster galaxies, the 'e(a)' galaxies (characterized by 
strong Balmer lines in absorption and emission lines of mod- 
erate intensity), are most easily interpreted as dust-hidden star- 
bursts. If the e(a) galaxies are the progenitors of the more com- 
mon k-Ha spectral-type galaxies, then it is likely that also k-i-a 
galaxies are dust-rich. 

Up to now, the available information on the IR properties 
of cluster galaxies is rather limited. Mid-IR (MIR, hereafter) 
fluxes have been published for members of the nearby clusters 



Virgo, Coma, and A 1367 (Bosel li et al. nWTllTWgl Quillen 
et al. 1999 Contursi et al. I2001> . as well as for the medium- 
distant clusters A370 (Soucail et al. 1999 Metcalfe et al. 2003 
hereafter Paper I), A 1689 (Fadda et al.;2000 Due et al. .2002j . 
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(Pierre et al. IT996b . A2218 (Barvainis et al. [19991 
Paper I), A2219 (Barvainis et al.[l999l, A2390 (Lemonon et 
al. 19981 Altieri et al. [T999l Paper I), and C10024-H1654 (Coia 
et al. 2004a C04 hereafter). These analyses have shown that, 
on average, optically-selected normal early-type galaxies have 
their MIR emission dominated by the Rayleigh-Jeans tail of 
the cold stellar component, while late-type galaxies have their 
MIR emission dominated by the thermal emission from dust 
(Boselli et al. .l998j . In the Coma cluster, the k-i-a galaxies were 
only found to have enhanced 12 yum emission when also op- 
tical emission-lines were present (Quillen et al. 1999). Since 
early-type galaxies are the dominant cluster population, one 
does not expect to find a large number of strong-IR emitters 
among cluster galaxies at wavelengths where the stellar pho- 
tospheric emission becomes negligible. Nevertheless, in C04 
we report the detection of a large number of very bright IR- 
galaxies in the distant cluster C10024-I-1654, an indication that 
the IR properties of cluster members can vary with redshift, and 
among different clusters. 

Outside clusters, there is now ample evidence for a signifi- 
cant evolution of the population of MIR-detected field galaxies 
(Altie ri et al. [79991 Aussel et al.[T999l Elbaz et a l.[T999l[2002[ 
lUUT; Paper I; Sato et al. HM^ Serjeant et al. HM). These 
MIR field galaxies have a redshift distribution peaked around 
z ~ 0.7. Their MIR-emission is mostly powered by starbursts 
(rather than AGNs), with typical star formation rates (SFRs) of 
~ 50 Mq yr"' . The integrated source counts of the MIR galaxy 
population exceed by an order of magnitude the predictions of 
no-evolution models based on the local MIR luminosity func- 
tion (Franceschini et al. 200 L Paper I). Models that explain 
the high MIR source counts require either a strong evolution 
of the whole luminosity function or a strong evolution in both 
luminosity and density of a sub-population of starburst galax- 
ies (Chary & Elbaz 2001 ; Franceschini et al UMU XuEOOO)- 
The integrated emission from the field MIR sources detected by 
ISOCAM accounts for most of the cosmic infrared background 
detected by Puget et al. ([T9961 see Elbaz et al. [2002b . 

In this paper we analyse the properties of IR-selected 
galaxies in the A2218 cluster, located at a mean redshift 
z = 0.175. We base our analysis on the data obtained in our 
gravitational-lensing deep-survey programme, which was part 
of the "Central Programme" of the ISO mission. These data 
were used in Paper I, in conjunction with ISOCAM data for 
other clusters, to derive deep source counts at 6.7 and 14.3 yum. 
The counts we derived in Paper I confirm and extend earlier 
findings of an excess by a factor of ten in the population of 
14.3 yum emitters with respect to no-evolution models (Altieri 
et al. 1999 Aussel et al. 1999; Elbaz et al. 1999), and they 
agree with the counts of 6.7 /zm sources reported by Sato et al. 
{"MT) and Aussel et al. C^^W). We complement the ISOCAM 
MIR data with photometric and spectroscopic data at optical 
and near-IR (NIR) wavelengths, taken from the literature. We 
fit the optical-MIR spectral energy distributions (SEDs) of the 
ISOCAM sources with the GRASIL models of Silva et al. 
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C1998' S98 hereafter) to constrain the spectral types, total IR 
luminosities, and SFRs of the cluster and background sources. 

In § 12 we describe our data-set. In § |5] we describe the 
model fits to the observed SEDs of the ISOCAM extragalac- 
tic sources, and in §|5we describe how we determine the IR 
luminosities and SFRs of these sources. In the two following 
sections we describe the results of our analysis, separately for 
cluster sources (§|5} and field sources (§|6j. We address the 
issue of the BO-effect in the IR in §0 In §11] we discuss our 
findings, and compare them with previous results from the liter- 
ature. In §|9]we summarize our results and provide our conclu- 
sions. Some aspects of the SED fitting procedure are detailed 
in AppendixlAl We discuss a few individual ISOCAM sources 
and their SEDs in AppendixiBl 

In this paper we use Hq - 70 km s ' Mpc"', Qq - 0.3 
and Qa = 0.7 (see, e.g., Melchiorri & Odman '2003V In this 
cosmology, the cluster luminosity distance is 845 Mpc, 1 ar- 
cmin corresponds to 178 kpc, and the age of the Universe at 
the cluster redshift is 11.3 Gyr. 



2. The Data 

The Abell 2218 field was observed with the LW2 and LW3 fil- 
ters of ISOCAM (Cesarksy et al. 1996 1, with reference wave- 
lengths of 6.7 and 14.3 //m, respectively, in the context of 
a gravitational-lensing deep survey programme, part of the 
"Central Programme" of the ISO mission. This programme is 
described in detail in Paper I. The total observing time was 6.2 
hours per filter. The observations were done in raster mode with 
the 3" per-pixel field-of-view, and a raster step size of 16", cov- 
ering a total area of 20.5 arcmin^, i.e. ^ 0.8 x 0.8 Mpc^ at the 
cluster distance. 

For a detailed description of the ISOCAM data reduction 
and source identification and reliability, we refer the reader to 
Paper I. The only modification with respect to Paper I concerns 
an improved photometry for the three extended sources 28, 30, 
and 45 (we use the ISOCAM source identification numbers of 
Paper I). As was pointed out in Paper I in the cases of sources 
28 and 45, the fluxes of these three sources as listed in that 
paper are underestimates, because their MIR photometry was 
derived under the point-source assumption, but these sources 
are extended. The total flux of these sources is between 1.5 
(source 45) and 2 (sources 28 and 30) times higher than quoted 
in Paper I. Since these three sources are cluster members, their 
flux corrections have no eff'ect on the field-source counts re- 
ported in Paper I. 

In total, 76 sources were detected in the A2218 field, with 
an average level of significance (i.e. the ratio of the source sig- 
nal to the 1-cr local noise floor) of 9 and 12 in the LW2 and 
the LW3 band, respectively. Four of these sources were also 
detected by Barvainis et al. ( 1999 i: their flux density estimates 
are in agreement with ours, within the errors. Among our 76 
MIR sources, 18 were detected in both bands, 30 in the LW2 
band only, and 28 in the LW3 band only. Source fluxes at 6.7 
and 14.3 yum are in the ranges 23-934 yuJy, and 90-919 juJy, re- 
spectively. Background sources have their fluxes amplified by 
the cluster gravitational lensing. After correction for lensing 



amplification (see Paper I), the faintest source detected at 6.7 
(14.3) ;um has an intrinsic flux of 5 (18)/iJy. 

We looked for optical and near-infrared (NIR) counterparts 
for the 76 ISOCAM sources, in the catalogues of Le Borgne 
et al. ([T9921 hereafter LPS), Small et al. (2001 hereafter SOI), 
and Ziegler et al. ( 120011 hereafter ZOl). The catalogue of LPS 
contains B, g, r, i, and z magnitudes for 729 objects in the cen- 
tral 4x4 arcmin^ of A2218. LPS also provide redshifts for 
66 objects in the same field. The catalogue of SOI contains 
photometry in the F450W, F606W, and F814W filters of the 
Hubble Space Telescope - HST herafter -, and in the K, fil- 
ter, for 81 galaxies in the central 2.5 x 2.5 arcmin^ of A2218. 
SOI also provide morphologies for most of their sources, and 
redshifts for 34 of them. Finally, the catalogue of ZOl con- 
tains U, B, V, I magnitudes and redshifts for 48 early-type clus- 
ter members in the central 2.5 x 2.5 arcmin^, as well as mor- 
phologies for 19 galaxies. 

We estimated the morphologies of another 17 optical coun- 
terparts from an HST image in the F8 14W band. The HST mo- 
saic of Abell 2218 was observed with WFPC2 in the F606W 
filter under the GO programme 7343 (PI: Squires). Most of 
the observations were conducted in March and May 1999, and 
the remaining in May 2000. It comprises 24 pointings ob- 
served in the CVZ (continue viewing zone) during 2 orbits. At 
each position 12 exposures of 700 sec were obtained totaUz- 
ing 8.4 ksec. The data were retrieved from the ST-ECF/ESO 
archive. The data reduction was done using standard IRAF 
and STSDAS scripts and was combined into one mosaic using 
Swarp (www.terapix.fr/swarp). More details on these data can 
be found in Hudelot et al. ("2004). For those ISOCAM sources 
that lie in the field covered by this HST image, we show in 
FiglDtheir MIR isocontours at 6.7 and 14.3 jim (in blue and 
red, respectively) overlaid on the optical image (sources with 
stellar counterparts - see below - are not shown). 

We search for optical/NIR counterparts within a circle of 
12" radius around each ISOCAM source, in order to match 
(roughly) the size of the ISOCAM PSF diameter at 14.3 //m 
(Okumura 1998 1. Following the method of Flores et al. ( I1999> . 
we determine quantitative likelihoods of all possible counter- 
parts. These likelihood estimates are based on the distance be- 
tween the ISOCAM source and its possible counterpart, and on 
the magnitude of the possible counterpart (see Flores et al."1999| 
for details). The mean (median) distance between the ISOCAM 
sources and their most likely counterparts is 1.4" (1.2"), and 
all sources, except two, have distances < 4" (sources 5 and 60 
have distances of ~ 7" from their most likely optical counter- 
parts). Visual inspection of the ISOCAM LW2 and LW3 maps 
overlaid on optical images (see Fig. ^ and Figs. 6 and 7 in 
Paper I) confirms the identification of the counterparts obtained 
by the maximum likelihood procedure. The likelihoods of the 
assigned optical/NIR counterparts of our ISOCAM sources are 
very high. Cases where confusion might play a relevant role are 
detailed in AppendixiBl 

Since the ISOCAM flux densities are total flux densities, 
we must convert the apparent magnitudes of the optical/NIR 
counterparts into total magnitudes. SOI provide total magni- 
tudes in the Ks-hand, and £450 - /^h, Vatd - hi4, hi4 - K, 
aperture colours, that we use to derive approximate total £450, 
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Fig. 1. MIR isocontours of ISOCAM sources, at 6.7 and 14.3 //m (in blue and red, respectively) overlaid on an optical HST image 
intheF814Wband. 



Vgoe, and /814 magnitudes (assuming negligible colour gradi- 
ents within the galaxies). Similarly, we determine approximate 
total U, B, and V magnitudes, using the /-band total mag- 
nitudes and U - I, B - I, V - I aperture colours in ZOl. 
We then convert ZOl's magnitudes into the photometric sys- 
tem of SOI, using the mean magnitude difference (in bands 
of similar central wavelength) of objects in common to the 
two data-sets. Finally, we transform LPS' isophotal magni- 
tudes into total, by using the following relations between LPS' 
and SOl's magnitudes, determined for objects in common to 

Blps - 0.60, If^l' = kps - L04, 
0.5 X (rips + gLPs ) - 0.50. Clearly, the magnitudes 
obtained using these empirical relations have larger errors than 
the original ones, because of the scatter in the relations (~ 0.2 
magnitudes), and these larger errors are taken into account in 
our analysis. 



the two data-sets: Bf^l 

\/total 

^ sm 



All the optical and NIR magnitudes are corrected for galac- 
tic extinction, and then converted into flux densities using the 
relations of Fukugita et al. ( 1995 1 and Tokunaga ( 2000). 

Nine of the 76 ISOCAM sources are identified with stars. 
These are sources no. 3, 4, 11, 12, 15b, 16, 56, 71, and 72. 
Seven of our ISOCAM sources are outside the surveyed areas 
of the LPS, SOI, or ZOl catalogues. These are sources no. 19, 
36a, 43, 48, 50, 55, and 73b. The lack of photometry in opti- 
cal/NIR bands precludes a thorough analysis of their properties. 
These 16 sources are not considered in this paper. We are thus 
left with a sample of 60 extragalactic ISOCAM sources with 
optical counterparts and available optical and/or NIR magni- 
tudes. 

3. The Spectral Energy Distributions 

We construct Spectral Energy Distributions (SEDs) for the 60 
non-stellar ISOCAM sources with available optical/NIR pho- 
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Fig. 1. Continued. 



tometry. When either the 6.7 or the 14.3 fj.m flux density is not 
available, an upper limit is used instead. This limit is equal to 
the 50% completeness limit for the region where the ISOCAM 
source is located (see Table 3 in Paper I), or twice that value 
for the extended sources 28, 30, and 45. With this definition, 
MIR upper limits approximately correspond to 4 cr limits. SED 
models that violate this limit are rejected in the fitting proce- 
dure. 

In order to fit the observed SEDs of our ISOCAM sources 
we consider a battery of models computed with the GRASIL 
code of S98. This code computes the spectral evolution of 
stellar systems by taking into account the eff'ects of dust. In 
particular, it takes into account several environments with dif- 
ferent dust properties and distributions, such as the AGB en- 
velopes, the diffuse interstellar medium, and the molecular 
clouds. GRASIL models consider diff'erent kinds of dust par- 
ticles, including the very small grains and the polycyclic aro- 
matic hydrocarbons that are mainly responsible for the MIR 



emission detected by ISOCAM. GRASIL models have been 
shown to provide excellent fits to the observed SEDs of nearby 
galaxies, in both quiescent and starburst phases, with infrared 



data obtained with ISO (S98; Mann et al. 120021 C04). FuH de- 
tails of the GRASIL model are given in S98. 

The available photometric data for our MIR-selected 
sources do not allow an accurate modelling of each individ- 
ual galaxy SED. Rather, we perform the comparison of the ob- 
served SEDs with a limited set of models meant to be represen- 



tative of four broad classes of spectral types. Mann et al. (I2002> 
have considered five model SEDs, while we have considered 20 
model SEDs in C04, i.e. 10 models each seen at two different 
ages. Here we consider 30 models, i.e. the same 10 models of 
C04, each seen at three different ages. Some of these models 
were taken directly from the GRASIL web site', some were 
kindly provided by L. Silva (private comm.), and some were 



' http://web.pd.astro.it/granato/grasil/grasiI.html 
or:_http://adlibitum.oat.ts.astro.it/silva/grasil/grasil.htnil] 
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Table 1. Results of the model SED fitting 



ISOCAM 
source id. 



Optical 
counterpart 



morphology 

(reference) 



Best-fit 
model 



Nsfr(O.l) 



Nsfr(l.l) 



Redshift 
(reference) 



Membership 



Quality 
of the fit 



ISO_A2218_02 


L664 


- 


Sa5 


0.4 


0.5 


0.53 (Paper I) 


f 


p 


IS0^22 18.05 


L621 


Sa 


Sa5 


0.4(0.1-0.8) 


0.1 (0.1-3.4) 


0.175 ±0.15 


c 


G 


ISO.A22 18.06 


L630 


E/SO 


Em5 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1802 (LPS) 


c 


G 


ISO.A22 18.08 


L587 


Sc 


Mps2 


1.4(1.1-3.0) 


1.3(1.3-1.4) 


0.68 (Paper I) 


f 


G 


ISO.A22 18.09 


L565 


SBa 


Sa2 


1.1 (0.4-1.8) 


1.2(0.5-2.4) 


0.68 (Paper I) 


f 


G 


ISO.A2218.10 


L557 


- 


Sb5 


1.1 (< 1.8) 


1.2 (< 3.4) 


0.00 ± 0.37 


- 


U 


ISOJ*l2218.13 


L537 


S 


M5 


3.0(1.1-3.0) 


1.3(1.3-1.4) 


0.90 ±0.16 


f 


G 


ISO.A2218.14 


L535/Z2076 


E 


Es2 


0.0 (< 0.3) 


0.00 (< 0.6) 


0.1827 (LPS) 


c 


G 


ISO.A2218.17 


L484 


Sa 


SblO 


0.6(0.4-1.1) 


0.6 (0.5-3.4) 


0.10 ±0.08 


c 


G 


ISO.A2218.18 


L482/S4013 


E/SO 


Em5 


0.0 


0.0 


0.1778 (LPS) 


c 


G 


ISOJV22 18.20 


L467/S159 


SO (SOI) 


Aps2 


0.5 (< 1.1) 


1.6(0.1-2.4) 


0.476 (E98) 


f 


G 


ISOJ*i22 18.21 


L436/S137/Z1516 


SBO/a (SOI) 


Em5 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1638 (SOI) 


c 


G 


ISO.A22 18.22 


L430/S337/Z2604 


E 


Em5 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1800 (LPS) 


c 


G 


ISO.A22 18.23 


Z1976 


E 


Em2 


0.0 (< 0.3) 


0.1 (<0.6) 


0.1686 (ZOl) 


c 


G 


ISO.A22 18.24 


L419/Z1142 


E 


Es2 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1641 (LPS) 


c 


G 


ISO.A22 18.25 


L409 


SBO/a 


Em2 


0.0 (< 0.3) 


0.1 (<0.6) 


0.1741 (LPS) 


c 


G 


ISO.A22 18.26 


L411/S103 


Sdm (SOI) 


Sa5 


0.4 (< 1.8) 


0.5 (< 3.4) 


0.175 ±0.20 


c 


G 


ISO.A22 18.27 


L395/S333 


Sc (SOI) 


Sa5 


0.4 (0.4-0.8) 


0.5 (0.5-3.4) 


0.1032 (LPS) 


f 


G 


ISO.A22 18.28 


L391/S301 


cD (SOI) 


EmlO 


0.0 


0.0 


0.1720 (LPS) 


c 


P 


ISO.A22 18.29 


L381 


SBO/a 


Sc2 


1.8(0.4-1.8) 


1.4 (0.5-2.4) 


0.521 (E98) 


f 


G 


ISO.A22 18.30 


L373/S280/Z1552 


SBO/a (SOI) 


SalO 


0.1 (<0.1) 


0.1 (<0.1) 


0.1776 (LPS) 


c 


G 


ISO.A22 18.32 


L348/Z2270 


Sa(ZOl) 


EslO 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1676 (ZOl) 


c 


G 


ISO.A22 18.33 


L341/S351/Z1662 


E(SOl) 


EmlO 


0.0 


0.0 


0.1637 (LPS) 


c 


G 


ISO.A22 18.34 


L323 


S 


SalO 


0.1 (< 0.3) 


0.1 (< 0.6) 


0.179 (E98) 


c 


G 


ISO.A22 18.35 


L317/S420 


Scd (SOI) 


Mps5 


1.1 


1.4 


0.474 (E98) 


f 


P 


ISO.A22 18.36b 


L296/Z1888 


Sab (ZOl) 


Em2 


0.0 (< 0.3) 


0.1 (<0.6) 


0.1760 (LPS) 


c 


G 


ISO.A22 18.38 


L289 


S/Irr 


Mps5 


1.1(1.1-3.0) 


1.4(1.3-1.4) 


1.033 (LPS) 


f 


G 


ISO.A22 18.39 


Z849 


E 


E15 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1652 (ZOl) 


c 


G 


ISO.A22 18.40 


L287 


Irr 


Sb2 


1.6 (0.4-3.0) 


1.4(0.5-2.4) 


0.702 (E98) 


f 


G 


ISOJ*i.22 18.42 


L275 


Sa 


Sb2 


1.6(1.1-1.8) 


1.4(1.2-1.5) 


0.45 (Paper I) 


f 


G 


ISO.A22 18.44 


L262 


Irr 


Sc2 


1.8 


1.4 


0.596 (E98) 


f 


G 


ISO.A22 18.45 


L244/S307/Z1437 


SO (SOI) 


EmlO 


0.0 


0.0 


0.1646 (LPS) 


c 


G 


ISO.A22 18.47 


L235/Z1175 


SBO/a (ZOl) 


EmlO 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1765 (LPS) 


c 


G 


ISOJ^22 18.52 


L206/S323 


Sc (SOI) 


Sb5 


1.1(0.4-1.1) 


1.2(0.5-1.4) 


0.55 (Paper I) 


f 


G 


ISOJV22 18.53 


L205/S368 


Sc (SOI) 


Sc5 


1.6 


1.5 


0.693 (E98) 


f 


P 


ISO.A22 18.54 


L196/S401/Z1466 


E(SOl) 


Em5 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1798 (LPS) 


c 


G 


ISO.A22 18.57 


L149/S1057 


E(SOl) 


EmlO 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1753 (LPS) 


c 


G 


ISO.A22 18.58 


L148/S428/Z1343 


SO (SOI) 


EmlO 


0.0 


0.0 


0.1830 (LPS) 


c 


G 


ISO.A22 18.59 


L145/S646 


Sdm (SOI) 


Mps5 


1.1(0.4-1.1) 


1.4(0.5-1.4) 


0.628 (E98) 


f 


G 


ISO.A22 18.60 


L116/S4010 


E(SOl) 


Sa5 


0.4 (< 1.1) 


0.5 (< 3.4) 


0.2913 (LPS) 


f 


G 


ISO.A2218.61a 


L119/S638 


SO (SOI) 


SalO 


0.1 (< 0.8) 


0.1 (< 2.4) 


0.175 ±0.12 


c 


G 


ISO.A2218.61b 


L113/S633/Z2702 


SO (SOI) 


EslO 


0.0 


0.0 


0.1738 (LPS) 


c 


G 


ISO.A22 18.62 


L118/Z1293 


E (ZOl) 


Em5 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1758 (ZOl) 


c 


G 


ISO.A22 18.66 


L77/S4007 


SBO(SOl) 


Es5 


0.0 (< 0.3) 


0.0 (< 0.6) 


0.1681 (LPS) 


c 


G 


ISO.A22 18.67 


L75/S4004 


- 


Es2 


0.0 


0.0 


1.5 


- 


P 


ISO.A22 18.68 


L62/S4009 


Sa(SOl) 


Sa5 


0.4 (< 1.1) 


0.5 (< 2.4) 


0.42 (Paper I) 


f 


G 


ISO.A22 18.70 


L58 


- 


SalO 


0.1 


0.1 


0.1703 (LPS) 


c 


G 


ISO.A22 18.73a 


L38/S4015 


Sc (SOI) 


MIO 


2.8 (< 11.9) 


1.0 (< 3.4) 


0.30 ± 0.22 


- 


U 



built by ourselves by running the GRASIL code (also freely 
available on the web). 

The models we consider are the following: 

E: three models of early-type galaxies, characterized by an ini- 
tial burst of star formation lasting 0.5, 1.0, or 2.0 Gyr, and 
by passive evolution thereafter. Depending on the duration 
of the initial burst, we label these models 'Es' (short), 'Em' 
(medium), and 'El' (long). These models have been used in 



Granato et al. ( 1200 It and reproduce the range of SEDs of 
ellipticals and SOs in the nearby universe. 
S: three models of disk galaxies (spirals), characterized by dif- 
ferent values of: i) the gas infall timescales, and ii) the ef- 
ficiency term in the Schmidt-type law (see § 2.1 in S98). 
The star formation histories (SFHs, herafter) of these mod- 
els gently increases with time up to a maximum, and then 
gently decreases. The maximum SFR occurs at ~ 1, ~ 3, 
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Age (Gyr) 
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Fig. 2. The star formation histories of the GRASIL models used 
in this paper. Time since galaxy formation is on the x-axis, and 
SFR (normalised to the maximum value) is on the y-axis. Note 
the logarithmic x- and y-scale in the top and bottom panel, re- 
spectively. Top panel: Es (solid line), Em (dash-dotted hne), 
and El (dashed line) models. Middle panel: Sa (solid Une), Sb 
(dash-dotted line), and Sc (dashed line) models. Bottom panel: 
Mps (solid line), and Aps (dash-dotted line) models. Models M 
and A are like models Mps and Aps, respectively, except that 
they are seen at the time of the starburst (the spike in the Figure) 
rather than 1 Gyr after it. See the text for a detailed description 
of all models. 



and ~ 9 Gyr after galaxy formation, for the three models 
here considered, which we label 'Sa', 'Sb', and 'Sc', re- 
spectively, because their SEDs reproduce those of nearby 
spirals of the Sa, Sb, and Sc types. 
SB: two starburst models that provide good fits to the observed 
SEDs of the moderate starburst galaxy M82, and of the 
strong starburst galaxy Arp220. We refer to these two star- 
burst models as 'M' and 'A', respectively. The starburst is 
characterized by an e-folding time of 0.05 Gyr, and involves 
~ 0.01 and ~ 0.1 of the total mass of the galaxy, for the M 
and A model, respectively. 
PSB: the two starburst models mentioned above ('M' and 'A'), 
but observed 1 Gyr after the starburst event; we label these 
'post-starburst' models 'Mps' and 'Aps'. 

The SB models we consider here are rather extreme ex- 
amples of the starburst phenomenon. We considered including 
NGC6090 as another example of starburst galaxy, of interme- 
diate IR luminosity between Arp220 and M82. However, after 
suitable normalisation, the SED of M82 is not different enough 
from that of NGC6090 to justify including an additional model 
in our analysis, given the relatively large error bars on our 
(rather faint) MIR fluxes. 
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Fig. 3. Top panel: The SFR averaged over the last 1.1 Gyr vs. 
the SFR averaged over the last 0. 1 Gyr, both normalised by the 
SFR averaged over the entire SFH, for the different GRASIL 
models, grouped in five classes. Middle panel: the distribution 
of the best-fit models for the ISOCAM cluster sources in the 
Nsfr(l.l) vs. Nsfr(O.l) plane. The shading in each box is pro- 
portional to the number of best-fit models. Bottom panel: same 
as the middle panel, but for ISOCAM field sources. 



The SFHs of the 10 models are shown in Figure |2] Three 
SEDs were built for each of these ten models, by stopping their 
star formation histories at ages of 2, 5, and 10 Gyr The SB and 
PSB model SEDs are obtained from the same models, stopping 
their star formation histories at the time of the starburst event, 
and, respectively, 1 Gyr later. 

Since these GRASIL models were produced in order to fit 
the observed SEDs of nearby galaxies, it is not clear what the 
physical meaning of, e.g., a 2 Gyr Sa model is. Such a model 
might, or might not, resemble a real Sa of that age, depend- 
ing on the way an Sa galaxy forms and evolves. It is therefore 
more meaningful to parametrize these models in terms of their 
current and recent SFR's. More specifically, we define the nor- 
malized star formation rates Nsfr(0. 1 ) and Nsfr( 1 . 1 ) as the SFR 
of the model, averaged over the last 0. 1 and, respectively, 1 . 1 
Gyr, divided by the SFR averaged over the entire galaxy SFH. 
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These two quantities measure the shape of the model SFH, 
rather than the intensity of the star formation process. Using 
Nsfr(O.l) and Nsfr(l.l) we can group our 30 models into five 
classes (see Figure|3 top panel). Values of Nsfr(l.l) < 1 char- 
acterize galaxies that have experienced their main star forma- 
tion activity more than 1 . 1 Gyr ago. Among these, those with 
Nsfr(O.l) > 0.1 have some residual star formation activity left. 
When Nsfr(l.l) > 1, we distinguish three cases, depending on 
the value of Nsfr(O.l). When Nsfr(O.l) is very high, the galaxy 
is undergoing a starburst. Intermediate values of Nsfr(O.l) char- 
acterize galaxies near a maximum of star forming activity in 
their SFH (see Figure|5J. Finally, when Nsfr(O.l) is compara- 
ble or even lower than Nsfr( 1.1), the galaxy might have recently 
suffered a starburst, but the starburst phase is now over These 
are the models that we call post-starbursts (PSB), although they 
still retain substantial amounts of star formation activity, as 
the starburst event has used up only a fraction of the available 
gas. Note that the range of Nsfr(O.l) and Nsfr(l.l) values that 
characterizes the PSB galaxies is also shared by other models, 
namely those of normal spirals observed shortly after the peak 
of their SFH. 

We compare the observed and model SEDs with a stan- 
dard ;t'^ procedure (see Appendix IaI for a detailed description 
of some technical aspects of our SED-fitting procedure). When 
the redshift of the optical counterpart to the ISOCAM source is 
known, there are two free parameters in the fit, the model SED 
and its normalization. When the redshift is unknown, it is taken 
as an additional free parameter We only accept the best-fit so- 
lution if the range of acceptable values for the source redshift 
is sufficiently narrow (see AppendixIXt. 

In order to constrain the fit to a galaxy SED and estimate 
its probability, the number of free parameters plus one must be 
less than the number of photometric data-points. This is true 
for 48 of the original 60 sources. For 40 of them, the redshift is 
known. 

Note that, currently, no AGN component is included in the 
GRASIL models. We do not expect to find many cases of AGN- 
dominated emission among MIR selected sources. Several 
studies (e.g., Franceschini et al. 2001 Rowan-Robinson 2001 ) 
have shown that most ISOCAM sources have their MIR emis- 
sion dominated by dust-reprocessed stellar radiation, especially 



Table 2. Additional ISOCAM sources 



at the lower flux levels. Elbaz et al. i2002> have found that only 
12 + 5% of the ISOCAM sources in the Hubble Deep Field 
North have their MIR flux dominated by dust re-processed 
AGN radiation. Among the 48 ISOCAM sources for which we 
try a SED-model fitting, we can therefore expect 6 + 2 in which 
the MIR emission is dominated by dust-reprocessed radiation 
from an obscured AGN. This number could be different if the 
relative number of AGNs among galaxies depends on the envi- 
ronment. 

The results of our SED-fitting analysis are summarized in 
Table [fl along with the optical counterparts and their mor- 
phologies when available. In the columns of Table ^ the fol- 
lowing information is listed: 

1 . ISOCAM source identification; 

2. identification of the optical/NIR counterpart in the LPS, 
SOI, and ZOl catalogues (L: LPS, S: SOI, Z: ZOl); 



ISOCAM 


Optical 


Redshift 


source id. 


counterpart 


(reference) 


IS0^22 18.07 


L599 


- 


ISO.A2218_15a 


L515 


- 


IS0^22 18.37 


L292 


- 


ISO.A22 18.41 


L283 


- 


ISO.A22 18.46 


L242 


0.654 (E98) 


ISO.A22 18.49 


L225 


- 


ISO.A22 18.51 


L220 


- 


ISO.A22 18.64 


LUl 


- 


ISO.A22 18.65 


L96 


0.64 (Paper I) 


ISO.A22 18.69 


L60 


^ 


ISO.A22 18.74 


L33 


- 


ISO.A22 18.75 


L26 


0.1688 (LPS) 



3. morphology of the optical counterpart, and reference (SOI, 
ZOl; if no reference is given, the morphology estimate is 
ours); 

4. best-fit model, followed by a number indicating the age of 
the model in Gyr; 

5. Nsfr(O.l), i.e. the SFR of the best-fit model, averaged over 
the last 0. 1 Gyr, divided by the SFR averaged over the en- 
tire SFH of the model. The 95 % confidence level (c.l. here- 
after) range (in brackets) is also given, when applicable, i.e. 
when other models provide an acceptable (j^'^ -probability 
of die fit > 0.05) fit to the observed SED, and when the 
Nsfr(0. 1 ) values of these models are different from the best- 
fit one; 

6. same as Col. 5, but for Nsfr( 1.1) (the SFR averaged over the 
last 1 . 1 Gyr, divided by the SFR averaged over the entire 
SFH),inlieuofNsfr(0.1); 

7. spectroscopic redshift, if available, and reference (LPS, 
SOI, ZOl, Paper I, and Ebbels et al. [D^ hereafter E98), 
or photometric redshift obtained through the SED fitting 
procedure, with its uncertainty (rms of the redshift values 
found among all the acceptable fits, see Appendix \K\. if 
available; 

8. cluster/field membership (labelled c/f, respectively), based 
on the spectroscopic or photometric redshift estimates of 
Col. 6; 

9. quality of the fit (G: good, ;t'2 -probability of the fit > 0.05; 
P: poor, ;^'^ -probability of the fit < 0.05; U: unconstrained, 
;f^ -probability of the fit > 0.05, but the photometric redshift 
solution is unconstrained). 

In Table 121 we fist the 12 ISOCAM sources without suffi- 
cient photometric data for a reliable SED-fitting analysis. In 
Col. 1 we list the ISOCAM source identification, in Col. 2 the 
identification of the optical/NIR counterpart in the LPS cata- 
logue, in Col. 3 the spectroscopic redshift, if available, and its 
reference. 

In summary (see also Table O, we obtain acceptable fits 
for 41 out of the 48 galaxy SEDs, five of which have no 
spectroscopic -z estimate. Of these five, four have photometric- 
z's consistent with the mean cluster z, so we consider them as 
cluster members. In total, acceptable fits are obtained for 27 
cluster members and 14 field galaxies. For two sources the best 
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Table 3. Summary of the SED fitting analysis 
Sample 



number 



Non-stellar sources 60 

Sources with well-sampled SED 48 

Sources without acceptable SED model fit 5 

Sources with unconstrained SED model fit 2 

Sources with acceptable SED model fit 41 

Cluster sources with acceptable SED model fit 27 

Field sources with acceptable SED model fit 14 



fit is unconstrained because of a degeneracy problem (different 
models at very different z's provide fits of comparable quality). 
For another five sources, no acceptable fit was found (they are 
discussed in Appendix|B}. 



4. The IR Luminosities and Star Formation Rates 

Computing the IR luminosities, L/r's, of our MIR sources re- 
quires knowledge of both their redshifts and SEDs. We can 
therefore determine the L/^'s for 41 ISOCAM sources in the 
A2218 field. 

We adopt two methods to determine Ljr: (a) We integrate 
the best fit model SED over the range 8-1000 jum; (b) we use 



the empirical relations of Elbaz et al. J2002I eqs. 13 and 14 in 
that paper), between Lir and the K-corrected monochromatic 
MIR luminosities at 6.7 and 14.3 fim. The best-fit model SED 
is used to estimate the K-correction. In both cases, the lumi- 
nosities of background sources are corrected for the lensing 
amplification factors (taken from Tables 5 and 8 of Paper I). 

The LjR estimates obtained by the two methods are very 
similar, except for those sources with SEDs that are best fit 
by E-type models, for which method (b) returns systematically 
higher luminosities than method (a) - see Fig. |4] This is not 
surprising, since the relations of Elbaz et al. J2002> are consis- 
tent with those of Chary & Elbaz ( 2001j, which were derived 
on a sample of galaxies with significant star formation activity. 



The SEDs of the galaxies in Chary & Elbaz's ( 1200 H sample are 
clearly different from those of E-type galaxies, where the MIR 
emission is dominated by the Rayleigh- Jeans tail of the cold 
stellar component (see, e.g., Boselli et al. ll998> . not by thermal 
emission from dust. Because of this problem, in the following 
we adopt the L/R-estimates obtained with method (a). We also 
calculate the NIR-luminosities, L^ir, using the same method, 
i.e. by integrating the best fit model SED over the range 0.9-3.0 
fj.m. 

We then compute the SFRs from Kennicutt's J 19981 1 rela- 
tion between a galaxy SFR and its Lir, 



SFR(Moyr-') 



1.71 X 10"'° L/r/Lq. 



(1) 



In Tablel^we list in Col. 1 the source identification (sources 
without spectroscopic redshift are listed in italic), in Cols. 2 
and 3, L^ir and Ljr, respectively, in units of 10'' Lq, and fi- 
nally in Col. 4 the SFR, in units of M© yr '. Only the 27 cluster 
sources and the 14 field sources with acceptable model fits to 
their SEDs are listed. 
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Fig. 4. The ISOCAM sources IR luminosities obtained from the 
best-fit model SEDs plotted vs. the IR luminosities obtained 
from the K-corrected LW3 luminosities, or, when these are not 
available, from the K-corrected LW2 luminosities, through the 
relations of Elbaz et al. C2002j . The dashed line is the identity 
line. Crosses indicate those sources whose SED is best fit by an 
E-type model. 



5. Cluster sources 

5. 1 . Internal properties 

In Fig. |5l we show the SEDs of the 27 cluster sources for 
which an acceptable fit was found (see Table|3}, with the best- 
fitting models overlayed. Most cluster sources are character- 
ized by a SED that decreases from the NIR to the MIR, typ- 
ical of passively evolving galaxies. Models with Nsfr(O.l) ^ 
Nsfr(l.l) ^ provide the best-fit to the SEDs of 20 cluster 
sources (see Fig.|3] middle panel). Models with slightly higher 
values of Nsfr(0. 1 ) and Nsfr( 1.1) (but still below unity) provide 
the best-fit to the SEDs of the remaining seven cluster sources 
(three of which however lack spectroscopic redshift determi- 
nations, so their cluster membership is not certain). Models 
with Nsfr(O.l) - Nsfr(l.l) ^ are still acceptable for four of 
these seven galaxies (see Table|3}, and only three of them have 
Nsfr(O.l) significantly greater than zero, implying significant, 
albeit small, ongoing star formation. 

In order to make a more efficient use of our data, we com- 
bine the 27 SEDs of the ISOCAM cluster sources. We assume 
the mean cluster redshift for all cluster sources, and we scale 
the observed SED of each source by the rest-frame //-band 
flux, as estimated from the best fit model. Not all sources are 
detected in the LW3 ISOCAM band. We estimate the aver- 
age LW3 flux using both the measured LW3 flux values and 
the LW3 flux upper limits. The average LW3 flux we derive 
is therefore an overestimate, and we treat it as an upper limit. 



10 
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Fig. 5. The observed SEDs of the 27 cluster sources for which we obtain an acceptable fit (dots with 2-cr error bars), with the 
best-fit GRASIL models overlayed (solid lines for sources with spectroscopic z, dashed line for sources with z determined from 
the model fitting procedure). The rest-frame wavelength, in fim, is plotted along the x-axis, and the flux density (in normalized 
units) is plotted along the y-axis. Upper limits are indicated by arrows. The approximate widths of the LW2 and LW3 bands are 
indicated. 



The resulting average SED is shown in Fig. |6l with the best 
fit model, an Em 10, superposed. This model (and all accept- 
able models, at > 5% c.l.) has Nsfr(O.l) = Nsfr(l.l) = 0. This 
analysis therefore confirms that, on average, ISOCAM cluster 
sources have negligible ongoing, or even recent, star forma- 
tion. This result indicates they must have formed the bulk of 
their stars quite some time before the observing epoch. 



Our conclusion is in agreement with what can be inferred 
from optical and NIR data. In fact, morphologically, most 
ISOCAM cluster sources are classified as ellipticals or SO's, 
which in general are characterized by old stellar populations. 
A similar conclusion can be reached from the analysis of their 
optical/NIR colours. InFig.0weplotthe V(,o6-h\4 vs. /g^-A'^ 
colour-colour diagram from the cluster galaxies in SOl's sam- 
ple. ISOCAM sources lie in the sequence defined by the red- 
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ISOCAM 


Lnir 


LlR 


SFR 


source id. 


lO'Le 


lO'Lo 


Moyr' 


Cluster sources 


ISOJl2218.05 


7.2 


17.3 


2.9 


IS0^22 18.06 


58.2 


0.7 


0.1 


ISO.A2218.14 


42.8 


0.8 


0.1 


ISO A22 18.17 


1.4 


4.1 


0.7 


ISO.A2218.18 


44.1 


0.6 


0.1 


ISO.A22 18.21 


18.8 


0.2 


0.0 


ISO.A22 18.22 


27.1 


0.3 


0.1 


ISO.A22 18.23 


32.3 


0.7 


0.1 


ISO.A2218.24 


31.2 


0.6 


0.1 


ISO.A22 18.25 


9.7 


0.2 


0.0 


ISOJi2218.26 


2.5 


5.9 


1.0 


ISO.A22 18.30 


55.0 


26.4 


4.5 


ISOJSl22 18.32 


15.6 


0.3 


0.0 


ISO.A2218.33 


40.2 


0.5 


0.1 


ISO.A2218.34 


6.4 


3.1 


0.5 


ISO.A22 18.36b 


27.0 


0.5 


0.1 


ISO.A22 18.39 


24.8 


0.4 


0.1 


ISO.A22 18.45 


83.4 


1.0 


0.2 


ISOJ^2218.47 


53.6 


0.7 


0.1 


ISO.A2218.54 


37.3 


0.5 


0.1 


ISO.A22 18.57 


50.7 


0.6 


0.1 


ISO.A2218.58 


42.7 


0.5 


0.1 


ISOJ^2218.61a 


16.5 


7.9 


1.3 


ISO.A2218.61b 


36.2 


0.7 


0.1 


ISO.A22 18.62 


33.6 


0.4 


0.1 


ISO.A22 18.66 


23.9 


0.4 


0.1 


ISO.A22 18.70 


9.8 


4.7 


0.8 


Field sources 


ISO.A22 18.08 


23.6 


353.3 


60.1 


ISO.A22 18.09 


25.5 


210.3 


35.7 


ISOJl2218J3 


41.6 


526.6 


89.5 


ISO.A2218.20 


6.1 


42.5 


7.2 


ISO.A22 18.27 


5.6 


13.4 


2.3 


ISO.A22 18.29 


7.9 


73.2 


12.4 


ISO.A22 18.38 


103.5 


736.5 


125.2 


ISO.A22 18.40 


2.4 


25.9 


4.4 


ISO.A2218.42 


11.8 


125.9 


21.4 


ISO.A2218.44 


4.5 


41.7 


7.1 


ISO.A22 18.52 


38.9 


225.6 


38.3 


ISO.A22 18.59 


22.8 


162.3 


27.6 


ISO.A22 18.60 


9.4 


22.5 
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der (and more luminous) galaxies, except source 26, a galaxy 
with an Sdm optical morphology (whose cluster membership, 
as mentioned before, is not certain). On the basis of this colour- 
colour diagram, SOI determine ages from 5 to 10 Gyr for 
the redder and more luminous A2218 cluster members (see 
Figure 3 in SOI). 

The consistent age estimates derived from the optical and 
the IR data suggest that there is little dust, on average, in the 
A2218 cluster galaxies, and that the MIR emission of ISOCAM 
cluster sources comes from the Rayleigh-Jeans tail of the pho- 
tospheric emission from cold stars (see, e.g., Boselli et al. 1998i) 
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1 10 

rest — frame X (/iin) 

Fig. 6. The average SED of the 27 ISOCAM cluster sources. 
Error-bars denote the rms among the flux values. The arrow in- 
dicates the upper limit on the LW3 flux value. The approximate 
widths of the LW2 and LW3 bands are indicated. The 27 SEDs 
were normalised to the //-band flux density, as estimated from 
the individual best fit models, before averaging them. The best 
fit EmlO model SED is also shown (solid line). 



The L/R- and SFR-distribution of cluster sources are shown 
in the two top panels of Fig.|8] The average L/s is 3 x lO'^L©, but 
the median is only 6 x 10** L©, and in fact most cluster sources 
have L/R < 10'' L©, corresponding to SFR < IMoyr '. There 
are a few sources with moderate L/r's and SFRs, but none qual- 
ify as a Luminous IR Galaxy (LIRG, L/r > lO^L©, see Genzel 
& Cesarskv l2000> . Source 30 has the highest SFR, but even in 
this case it is only ~ AMq yr"'. 

If the galaxy is small, even a relatively small SFR can be 
considered a sign of strong ongoing star formation activity. In 
order to investigate this, we consider the LirJLmir ratio vs. 
Lnir- The SFR scales linearly with L/r (eq. [0, and Lmir is 
a fair indicator of the total baryonic mass of a galaxy (e.g. 
Gavazzi et al. ll996> . so LirILnir is proportional to the SFR 
normalised by the baryonic mass of a galaxy. We plot LirIL^ir 
vs. Lnir for cluster sources in Fig. |9] (top panel). For most 
cluster sources Ljr/Lnjr < 0.03, and for all spectroscopically- 
confirmed cluster members Ljr/Lnir < 0.5. Only three non- 
spectroscopically confirmed cluster members (sources 5, 17, 
and 26) have Ljr > L^m, and for one of them (source 5) the 
assigned optical counterpart is not certain (see § |2j. Thus, in 
general, even the normalised SFRs of ISOCAM cluster sources 
are small, although the SFR per unit baryonic mass seems to be 
higher for galaxies of lower NIR luminosity. 

We conclude that most ISOCAM cluster sources are pas- 
sively evolving galaxies. Their MIR emission is of stellar 
photospheric origin. There are nonetheless ISOCAM cluster 
sources with non-zero, albeit small, SFRs, and the most active 
of these sources are those with the smallest baryonic masses. 
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Fig. 7. The Veoe - /si4 vs. /gn - Ks colour-colour diagram for 
cluster sources with data from SOI. ISOCAM galaxies are plot- 
ted as diamonds. 



5.2. Spatial and velocity distributions 

In order to obtain useful information from the spatial and ve- 
locity distributions of ISOCAM cluster members, we compare 
them with those of optically selected cluster sources. In Fig. 1101 
we show a smoothed projected number-density map, (obtained 
with the method of the adaptive kernel, see, e.g., Biviano et 
al. .1996 1, of the 266 likely cluster members from the LPS 
catalogue. The membership has been established by selecting 
galaxies within +0.25 magnitude of the colour-magnitude se- 
quence in the B-r vs. r diagram (see, e.g., Biviano et al. 1996 1. 
Superposed on the density map are the ISOCAM galaxies be- 
longing to the cluster, and with photometric data from LPS. 
The ISOCAM sources seem to share the same spatial projected 
distribution of optically-selected sources, except perhaps that 
they seem to avoid the main density peak. However, the relative 
fractions of optically-selected and MIR-selected cluster mem- 
bers contained within the highest-density cluster region (i.e. 
within the second highest isocontour of Fig.llO> are not signif- 
icantly different (23/266 and 2/27, respectively). And in fact, 
the whole spatial distributions of the two samples are not sig- 
nificantly different (22 % c.l.) according to the bi-dimensional 
Kolmogorov-Smirnov test (see, e.g., Fasano & Franceschini 
[T9871 . 

We now check the relative distributions of optically- 
selected and MIR-selected cluster members in velocity space. 
We select 79 cluster members from the full sample of 101 
galaxies with measured redshifts in the A2218 field, using the 
method of den Hartog & Katgert ( 1996 ). The velocity distribu- 
tion of the A2218 cluster members is shown in Fig. II II (note 
that velocities are rescaled to the cluster rest-frame using the 
average cluster velocity). The dashed and solid histograms re- 
fer to the 79 cluster members, and, respectively, to the sub- 
sample of 25 ISO cluster members with available redshift (see 
Tables [2 and |2ji. We run the Kolmogorov-Smirnov test to test 
the null hypothesis that the two samples have the same par- 
ent velocity distribution; the null hypothesis cannot be rejected 
(17 % C.L). 



The average velocity and velocity dispersion of the 79 
optically-selected cluster members (corrected for cosmologi- 
cal effects and velocity measurement errors, see Danese et al. 
[T980I . ai-e 52376 + 188 km s"' and 1412+[i^ km s"', respec- 
tively (here and hereafter, the biweigth estimators for the loca- 
tion and scale of a distribution are used, see Beers et al. ll990> . 
These values are in agreement with those given by Cannon et 
al. J1999> and Girardi & Mezzetti (2001 1, and are not signifi- 
cantly different from the average velocity and velocity disper- 
sion of the 25 ISOCAM cluster members (52091 + 388 km s^' 
and 1606^2^5 km s"', respectively). 

We conclude that the projected phase-space distribution of 
ISOCAM-selected cluster members is not significantly differ- 
ent from that of the global cluster population. 

It is nonetheless interesting to analyse the distribution of 
cluster galaxies in some more detail, to search for the existence 
of substructures. We find that the velocity distribution of the 79 
central cluster members is marginally different from a Gaussian 
(98 % C.L), according to the Anderson-Darling test (see, e.g., 
D'Agostino ll980> . even if a battery of other tests (skewness, 
kurtosis, tail-index, see Bird & Beers 1993 1 fail to find a signifi- 
cant difference. We also note that the velocity of the cD galaxy 
is significantly offset from the cluster mean (as indicated by 
the test of Gebhardt & Beers 1991 see also Fig.llH. The non- 
Gaussian velocity distribution of the cluster members, and the 
fact that the cluster cD is not at rest at the bottom of the clus- 
ter potential, suggest that the cluster is not fully dynamically 
relaxed (see, e.g., Girardi & Biviano .2003j . 

We then look for the presence of substructures in the 3d- 
space of positions and velocities. We do not detect significant 
substructure with the classical test of Dressier & Shectman 
(I1988> . However, a detailed modelling of the gravitational- 
lensing properties of the cluster has suggested the presence 
of two main mass concentrations, one centered on the cD and 
the other on the second brightest member, LPS244 (Kneib et 
al. 19951. The model for the mass distribution of Kneib et al. 
( 1995j is remarkably similar to the projected galaxy distribu- 
tion shown in Fig.[TOl(compare with Fig. 5 in Kneib et aL ll995> . 

Since there is evidence for a bimodal cluster structure, one 
could ask whether the ISOCAM cluster members belong to one 
structure rather than to the cluster as a whole. However, neither 
from the projected spatial distribution (see Fig. I10> . nor from 
their velocity distribution (see Fig. Ill> . is there any evidence 
for significant subclustering of the ISOCAM cluster popula- 
tion. 



6. Field sources 

In Fig. ^] we show the SEDs of the 14 field sources for which 
an acceptable fit was found (see Table |3}, with the best-fitting 
models overlayed. The SEDs of field sources are quite differ- 
ent from those of cluster sources, and are charachterized by an 
increasing flux density with wavelength. Such a SED is char- 
acteristic of actively star-forming galaxies, where part, or even 
most, of the stellar radiation is reprocessed by dust and re- 
emitted in the IR. In fact, models that provide the best-fits to 
the field-source SEDs have median values of Nsfr(O.l) - 1.1 
and Nsfr(l.l) = 1.4, indicating substantial ongoing and recent 
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Fig. 8. The distributions of L/r and SFR for ISOCAM cluster and field sources. Solid (dashed) histograms show the distributions 
for sources with (respectively, without) known redshift. Top left panel: the distribution of L/r for cluster sources. Bottom left 
panel: the distribution of Ljr for field sources. Top right panel: the distribution of SFR for cluster sources. Bottom right panel: 
the distribution of SFR for field sources. 



star formation, at a rate higher, on average, than the mean SFR 
over the galaxy SFH. Taking into consideration all the accept- 
able models, we note that Nsfr(O.l) and Nsfr(l.l) are both sig- 
nificantly higher than zero in 1 1 (out of 14) field sources. 

The ongoing star formation activity in field sources is con- 
sistent with them being mostly (late) spirals or irregulars, as 
expected from the morphology-density relation (e.g. Dressier 

The LjR of field sources correlates with redshift (see 
Fig. I13>. as expected in a flux-limited survey. In Fig. 1131 we 
show that the lower envelope of the observed Lju vs. z rela- 
tion approximately corresponds to the expected relation for a 
source with an observed 14.3 fj.m flux density 0.121 mJy (the 
50% completeness Hmit of the A2218 ISOCAM survey, see 
Paper I), corrected by the average amplification factor (0.6) due 
to the cluster lensing, and using the average K-corrections for 
the models that best fit the observed field source SEDs. 

The LjR- and SFR-distribution of field sources are shown in 
Fig. |8l(bottom panels). The average Ljr is 1.9x 1O"L0, and the 
median is 1.3 x IO^'Lq. Eight of the 14 field sources classify 
as LIRGs (~ 60% of all field sources). All LIRGs are sources 



with known redshift. These high L/jj's translate (via eq.[3 into 
high SFRs, 2-125 Moyr"', with an average (median) of 33 
(22) Moyr"'. The highest-SFR source in our sample is no. 38, 
a lensed object at z = 1.033 already detected in the MIR by 
Barvainis et al. J1999> . 

The distribution of the field sources in the Nsfr(l.l) vs. 
Nsfr(O.l) diagram (see Figure|3| bottom panel) can help us un- 
derstand better the nature of these field sources. Most of them 
are located in the central part of the diagram, characterized 
by similar values of Nsfr(l.l) and Nsfr(O.l), both above unity 
(see also Table [l), and only for one source is Nsfr(O.l) sig- 
nificantly higher than Nsfr(l.l). This suggests that we are not 
observing these galaxies during an exceptional starburst event, 
and that their current SFR was maintained at a similar level 
for quite some time. If anything, the median value of Nsfr( 1.1) 
of these 14 field galaxies is higher than their median value of 
Nsfr(O.l), implying a higher SFR in the recent past than to- 
day. Consistently, almost half of the galaxies in the field sam- 
ple are located in the Nsfr(0.1)-Nsfr(l.l) part of the diagram 
also occupied by PSB models (see Figure |3j bottom panel). 
However, on a galaxy per galaxy basis, the values of Nsfr(O.l) 
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Fig. 10. The projected space distribution of ISOCAM cluster 
sources (filled squares: spectroscopically confirmed members; 
open squares: sources whose cluster membership is established 
from the best fit model SED), and of optically-selected cluster 
sources (isodensity contours). The two crosses indicate the po- 
sition of the cD (near the centre) and of LPS244, the second 
brightest cluster galaxy. 
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Fig. 9. The ratio of the IR luminosity to the NIR luminosity, 
Lir/Lnir as a function of Lnir, for ISOCAM sources. Open 
symbols indicate sources without spectroscopic redshift esti- 
mates. Top panel: cluster members. Sources with Nsfr(O.l) > 
are crossed. Bottom panel: field galaxies. 



and Nsfr(l.l) are not significantly different, and other models, 
implying no major starbursts (e.g. intermediate-age spirals) fit 
the field galaxy SEDs as well as the PSB models. 

With the present data we are therefore unable to say 
whether the IR-selected field galaxies in our sample experi- 
enced one (or several) starburst events in the recent past, or 
whether they underwent a smoother evolution. However, the 
SED fitting analysis suggests that their SFRs have not changed 
significantly in the last ~ 1 Gyr. This conclusion implies that, 
on average, ~ 3 x 10"^ M© of stars were produced in the last Gyr 
before the observing epoch. This represents a substantial frac- 
tion of the total baryonic mass of these galaxies (as estimated 
from the best-fit GRASIL models) ~ 30%, on average. 

These ISOCAM field sources are therefore powerful star- 
forming galaxies, and they shine in the IR because a substan- 
tial fraction of their stellar luminosity is re-emitted in the IR. 
The amount of extinction is ~ 0.5-1.5 and ~ 1-2 magnitudes, 
respectively, in the rest-frame B and U bands, for the typical 
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Fig. 11. The rest-frame velocity distribution of optically- 
selected cluster members (dashed line) and of ISOCAM cluster 
sources (solid line). Arrows indicate the location in the veloc- 
ity space of the cD galaxy (at negative velocity) and the 2nd 
brightest member, LPS244. 



best-fit models. All field sources have L[r > IL^jr, and a few 
have LiRJLjsiiR > 10 (see Fig.|9l bottom panel). 

7. The Butcher-Oemler effect 

In Fig. [2] we show the B - r \s. r colour-magnitude dia- 
gram for all the sources in the LPS catalogue. Among these, 
the sources detected by ISOCAM are plotted as filled sym- 
bols, and confirmed cluster sources are plotted as diamonds. 
The dashed line is a best-fit to the main sequence of cluster 
members. All ISOCAM cluster sources lie close to the main 
colour-magnitude sequence. 

The fraction of blue BO galaxies was defined by Butcher 
& Oemler (1984) as the fraction of galaxies 0.2 mag bluer 
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Fig. 12. The observed SEDs of the 14 field sources for which we obtain an acceptable fit (dots with 2-cr error bars), with the 
best-fit GRASIL models overlayed (solid lines for sources with spectroscopic z, dashed line for the source with z determined 
from the model fitting procedure). The rest-frame wavelength, in fim, is plotted along the x-axis, and the flux density (in arbitrary 
units) is plotted along the y-axis. Upper Umits are indicated by arrows. The approximate widths of the LW2 and LW3 bands are 
indicated. 



than the early-type galaxies, and brighter than a K-corrected 
(but not evolutionary corrected) My = -20 (they used Hq = 
50km s"' Mpc"', and qo - 0.1). This absolute magnitude limit 
corresponds to an apparent magnitude limit r - 20.27, after 
correcting for the Galactic absorption in the direction of A2218 
(as given by Schlegel et al. 1998 1, and for the K-correction (we 
take the values from Poggianti .l997j . and using the V - r color 
of an early-type galaxy at z ~ 0.175 (Fukugita et al. ll995> . 



At the cluster redshift, a B - V colour difference of 0.2 mag 
corresponds roughly to the difference between the colours of 
an E and a SO galaxy; at the same redshift, the B - r colours 
of the same galaxy types differ by ^ 0.3 mag (Fukugita et al. 



I1995I I. Therefore, in order to compute the BO effect in A2218, 
we consider the galaxies with r < 20.27, and, among these, 
those bluer by > 0.3 mag than the colour-magnitude sequence 
in the B - r\s. r colour-magnitude diagram. 
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Fig. 13. The total IR luminosity as a function of redshift for 
ISOCAM field sources with spectroscopic redshift estimates, 
detected in the LW3 band. The solid line is the luminosity cor- 
responding to a flux of 0.121 mJy (the 50% completeness limit 
of the A2218 ISOCAM survey), reduced by a factor 0.6 (the 
average lensing correction factor), and K-corrected using the 
average correction for the models that best fit the observed field 
galaxy SEDs. 

Using the LPS data we estimate a BO galaxy fraction 
fs - 0.13 + 0.04, intermediate between the value given by 
Butcher & Oemler ([T984t . fg = 0.09 + 0.03, and that given by 
Rakos et al. (200l|l, fg = 0.23. However, if we only consider 
the confirmed cluster members, the fraction of blue BO galax- 
ies drastically decreases to fs - 0.07 + 0.04. Similar fractions 
of blue galaxies are found when selecting ISOCAM galaxies. 
When we consider all ISOCAM galaxies, the blue fraction is 
/b = 0.16 + 0.08; only five ISOCAM galaxies contribute to the 
blue fraction in the A2218 field (sources 2, 15a, 27, 61a, and 
64). Of these five galaxies, two are in the field, and three are at 
unknown redshift, although one, sources 61a, is a likely cluster 
member, based on the results of the SED fit analysis. 

In A2218 the BO effect is the same when galaxies are se- 
lected in the optical or in the MIR. Hence we would conclude 
that there is no MIR-BO eff'ect in A2218. However, we might 
have underestimated the size of the effect because of the rela- 
tively small field covered by ISOCAM observations. Ellingson 
et al. (I2OOI1I have in fact shown that in order to detect the BO 
effect it is necessary to sample radii larger than 0.5 r2oo, which 
is much beyond the radius of our ISOCAM observations of 
A2218. 

8. Discussion 

8.1. Cluster sources 

We have fitted the SEDs of 27 ISOCAM cluster sources (see 
Table O with GRASIL models. Of these 27 sources, most are 
characterized by a SED that declines from the NIR to the MIR, 
which is a typical feature of quiescent galaxies with little or no 
star formation. The SEDs of 20 ISOCAM cluster sources are 
best fit by models without episodes of star formation in the last 
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Fig. 14. The B - r vs. r colour-magnitude diagram for galax- 
ies in the A2218 field (magnitude and colour data are from 
LPS). Filled symbols represent ISOCAM sources. Diamonds 
represent cluster sources (big diamonds: spectroscopically con- 
firmed cluster members; small diamonds: cluster membership 
established through the SED-fitting procedure). X's identify the 
ISOCAM sources with a SFR > 4Mq yr '. The dashed line is 
the best-fit cluster colour-magnitude relation. Solid lines de- 
limit the region of the diagram where blue BO galaxies are lo- 
cated. 



1 . 1 Gyr For simplicity, hereafter we refer to these 20 galaxies 
as the passive cluster sample. 

The MIR emission of these galaxies corresponds to the 
Rayleigh-Jeans tail of the photospheric emission from cold 
stars, as found by Boselli et al. (1998 1 for early-type galax- 
ies in the Virgo cluster. The MIR flux densities of the galaxies 
detected by Boselli et al. (2003 1 in Virgo span a wide range 
(~ 0.1-1000 mJy); the 50% completeness limits of our survey 
of A2218 at 6.7 yum (see Paper I) corresponds to 130 mJy at 
the distance of Virgo, hence we have only detected the bright 
end of the distribution of early-type galaxies in A2218 (see also 
Figs.l7landll4>. In general, the optical morphologies of the pas- 
sive cluster members show they are early-type galaxies. 

Our ISOCAM-selected passive galaxies seem to share the 
properties of the average A2218 population of bright galaxies, 
for which SOI estimated an age of 5-10 Gyr. An age of 5 Gyr 
implies that the last major episode of star formation in these 
galaxies occurred at z ^ 0.9. 

The remaining seven ISOCAM-selected cluster sources 
show some evidence of a (small) star formation activity. We re- 
fer to these galaxies as the active cluster sample. Their SFRs 
are not strong, < 5 Mq yr"', and, according to the best-fit 
GRASIL models, are significantly lower than the average SFRs 
over the whole SFHs of these galaxies. However, galaxies of 
the active sample are among the faintest ISOCAM-detected 
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members of the A2218 cluster in the NIR (see Fig.|5] top panel; 
the active cluster sources are indicated with crosses). The ac- 
tive cluster sources have IR luminosities that are comparable 
to, or even higher than, their luminosities in the NIR. Our re- 
sults are therefore consistent with the general trend found by 
Boselli et al. (2003) in Virgo, namely that less massive cluster 
galaxies have higher ^14.3/^6.7 flux density ratios. 

Morphologically, most of the active cluster galaxies are not 
early-type, so it is perhaps not surprising that we find residual 
traces of recent star formation activity in their SEDs. However, 
two of them are SO-SBO/a. This is also not surprising, given the 
results obtained by SOI, who showed that 30% of the SO galax- 
ies in A2218 have relatively blue colours, indicative of a recent 
(or ongoing) star formation activity. This supports a scenario in 
which part of cluster SOs are formed via transformation from 
field spirals that have recently entered the cluster environment 
(see, e.g., Poggianti 2003 1. We must caution the reader that four 
galaxies of the active cluster sample are not spectroscopically 
confirmed members, and it is therefore possible that some of 
them are in fact field galaxies that we erroneously assign to the 
cluster (for one of them, source 5, even the assigned optical 
counterpart is not certain, see §|2j. 

The spatial and velocity distributions of the MIR cluster 
sources are not different from those of the overall A22 1 8 galaxy 
population. There is no evidence that they belong to substruc- 
tures of A2218 or to an infalling population. Rather, it seems 
that the MIR selection of A2218 cluster sources is not biased. 
Probably, this happens because our MIR cluster sample is dom- 
inated by the passive population selected at 6.7 fixn, where the 
emission is of stellar photospheric origin. We would have ob- 
tained a different view of the cluster MIR population had our 
6.7 fim observations been less deep. If we only consider the 
nine MIR cluster members detected at 14.3 jum, we find that 
only two of them belong to the passive population. It is there- 
fore important, when comparing A2218 with other clusters, to 
account for the surveys different sensitivities and different ar- 
eas. 

In this respect, it is particularly interesting to compare the 
number of cluster LIRGs detected by ISOCAM in five differ- 
ent clusters, namely A1689 (Fadda et al. 2000, Due et al. 2002), 
C10024+1654 (C04), A370 and A2390 (Paper I), and A2218. 
The advantage of using LIRGs to compare the number of IR 
active galaxies in different clusters is that a LIRG luminosity 
is high enough to lie above the limiting sensitivities of the sur- 
veys of all these clusters. In C04 we have recently determined 
the number of LIRGs for these five clusters, normalised by the 
surveyed areas and cluster masses. The fraction of LIRGs in 
C10024-1- 1654 turns out to be significantly higher than the frac- 
tion in A370, and than the average fraction in A1689, A2218, 
and A2390. There is however no significant difference in the 
LIRG fractions of the latter three clusters. 

In this same context it is also instructive to analyse the BO- 
effect in A2218, compared to the BO-effect in C10024+1654 
and A1689. In C10024-H1654 two thirds of the MIR cluster 
sources are part of the blue BO galaxies, and the MIR clus- 
ter sources contribute a third of all BO galaxies. In A1689, all 
BO galaxies are detected in the MIR (see Fig. 7 in Due et al. 



I2OO2I 1. In A2218, only one MIR cluster source qualifies as a 
blue BO galaxy, out of four cluster members. 

The reason for the different LIRG fractions and MIR BO- 
effects among different clusters is not clear. Part (but not all) 
of the differences could be due to the different surveyed ar- 
eas. While the LIRG fractions were computed after normalis- 
ing for the different surveyed areas, such normalisation might 
not be sufficient if the LIRG fraction significantly increases 
with the distance from the cluster centre. Such a dependence 
is suggested by the analysis of C10024-I-1654 (see C04), but 
more data are needed to confirm it. On the other hand, the ra- 
dial dependence of the BO-effect is an obvious consequence 
of the morphology-density relation (e.g. Dressier 1980 1, and 
larger physical areas have been covered by ISOCAM observa- 
tions in C10024+1654 and A 1689, than in A2218. Certainly, 
a larger physical area of A2218 needs to be covered in the IR, 
before we can make a soUd statement about the (lack of a) MIR 
BO-effect. After all, even in the optical the BO-effect remains 
undetected when only the cluster central regions are observed 



(Ellingson et al. l2001t . as in the present study. However, given 
that all BO galaxies are MIR emitters in A 1689, the fraction 
of MIR emitters among BO galaxies would not change if the 
surveyed area was smaller, and so it would still remain sub- 
stantially higher than the corresponding fraction in A2218. 

The different LIRG fraction is not simply related to the 
average evolution of clusters with redshift, since A370 and 
C10024-I-1654 have similar redshifts and yet have very dif- 
ferent LIRG fractions. In C04 we suggest that an ongoing 
merger could be at the origin of the starbursting activity 
in the C10024-H1654 galaxies. However, also A2218 is un- 



dergoing a merger with a substructure (Girardi et al. 119971 
Cannon et al. 1999 Machacek et al. 2002 ). An ongoing cluster- 
subcluster merger can stimulate starbursts in cluster galaxies 
through the tidal effects of a rapidly changing gravitational 
field (see, e.g., Bekki 2001 ), but only if these galaxies are gas- 
rich. Differences in the average cluster populations before the 
merger occurs could explain the different impact of a cluster- 
subcluster merger on the stimulation of star-formation and of 
IR-emission in cluster sources. As a matter of fact, the fraction 
of emission-line galaxies in C10024-I-1654 (40 % according to 
Czoske et al. l2001> . is much higher than in A2218, and in rich 
clusters in general (Biviano et al. '1997^, and, correspondingly, 
the fraction of blue galaxies is quite low in A2218, possibly 
even lower than the average fraction in nearby clusters (Rakos 
et al. l200lll . 

We might thus speculate that the enhancement in the frac- 
tion of MIR active cluster members is not simply the conse- 
quence of an ongoing cluster-subcluster collision, but also of 
the integrated history of the cluster, from which it depends how 
many cluster galaxies still retain significant amounts of gas to 
fuel the star formation activity. However, a larger physical area 
coverage of galaxy clusters in the IR is needed before we can 
make a conclusive statement about the variation (if any) of the 
number of MIR galaxies among different clusters. 
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8.2. Field sources 

We have fit the SEDs of 14 field sources(see Table ^ with 
GRASIL models. These SEDs are very different from those 
of cluster sources, and are characterized by an increasing flux 
density with increasing wavelength from the UV-optical to the 
MIR. Models of star-forming galaxies without an AGN compo- 
nent provide satisfactory fits to the observed SEDs. Although 
we cannot exclude that an AGN component is also present, it 
is unlikely to be dominant for these 14 sources. The SEDs of 
three field sources (see Appendix iBt could not be fit with any 
of the models we considered. Possibly this is due to a problem 
of source confusion, or, possibly, the emission of these three 
sources is AGN-dominated. In the latter case, the fraction of 
MIR sources in our sample with an AGN-dominated emission 
would be ~ 0.2, comparable to the fractions reported in the lit- 



erature (Fadda et al. 120021 Elbaz et al. 20021. Characterizing 
the MIR SED of AGNs is however beyond the scope of this 
paper. 

Models providing acceptable fits to the 14 field galaxy 
SEDs have significant recent and ongoing star formation. To be 
specific, we find that the median SFR of these 14 field galaxies 
is 22 Moyr ', and that this corresponds to 1.1 the SFR aver- 
aged over the galaxies entire SFHs, according to the best-fit 
GRASIL models. According to the same models, the SFRs of 
tiiese galaxies have not changed much over the last ~ 1 Gyr, 
so that ~ 30% of the galaxies total baryonic mass was con- 
verted in stars in the last ~ 1 Gyr before the observing epoch 
(i.e. since z ~ 0.8, taking the median redshift of these sources, 
z - 0.6, as representative of the sample). Such an intense star 
forming activity can occur in a continuous smooth fashion or 
in a single starburst. In the former case, we would identify our 
MIR-selected field galaxies with massive (~ 10" Mq) spirals 
near the peak of activity in their SFH. In the latter case, the 
starburst is likely to have occurred some time before observa- 
tion, as SB models do not in general provide best-fits to the 
observed field galaxy SEDs, and, in fact, PSB models are pre- 
ferred (see Tableland Figure|31 bottom panel). Note however 
that since we only considered two SB and two PSB models, 
the evidence in favour of PSB vs. SB models is only tentative, 
but a full statistical approach of the kind described by Chary & 
Elbaz ("2001^ is not justified because of the rather large errors 
on our (quite faint) MIR fluxes. 

Among the two rather extreme SFHs described above, a 
more realistic one could be that described by Elbaz et al. 



Il2004> . in which the galaxies undergo a series of Luminous 
IR Phases (LIRPs). Such LIRPs would be produced by minor 
starburst events triggered by galaxy-galaxy interactions, but not 
necessarily major collisions or major mergers. This scenario is 
supported by the observations of Flores et al. J1999l i who have 
found that the SEDs of many ISOCAM galaxies selected in 
the Canada-France Redshift Survey are not indicative of young 
dusty starbursts, but, rather, of normal spirals with a significant 
population of A stars, suggesting previous (recent) starburst- 
ing episodes. Similarly, Mann et al. ( 2002 1 have shown that the 
SEDs of several ISOCAM sources in the Hubble Deep Field 
South are adequately fitted by models of normal spirals, rather 
than by starburst models. 



The morphologies of the field ISOCAM galaxies in our 
sample span a wide range, from E to Sdm and Irr, but late- 
types dominate. Peculiarities and asymmetries are visible in a 
few cases, possibly indicating past or ongoing interactions with 
other galaxies, which could be partly responsible for an en- 
hanced star formation activity (Condon et al. I1982> . Previous 
studies have already suggested that interactions might be quite 
common among ISOCAM-selected field galaxies (Flores et al. 
[T9991 Elbaz & Cesarsky 2003 Hammer et al. l2004l . 

The LIRG fraction among our MIR field sources (~ 60%), 
and their average SFR (33 Mq yr"') are in good agreement with 
the values found in other ISOCAM field surveys (e.g. Chary & 
Elba z ,200L , Elbaz et al. 2002. Elbaz & Cesarski .2003. Mann 
et al.l2002>. In this context it is useful to warn the reader about 
the different nomenclature used here, with respect to the one 
adopted by Chary & Elbaz (2001 1 and Elbaz et al. 



oocm . 

We all use the terms LIRGs and ULIRGs to indicate, respec- 
tively, galaxies with 10" < L/r/L© < 10'^ and L/r/L© > lO'^. 
However, Chary & Elbaz (2001.) and Elbaz et al. QQ02 1 use the 
term 'starburst' to indicate galaxies with intermediate IR lumi- 
nosity (10'° < Ljk/Lq < lO"), while we use the same term 
to indicate that the galaxy is undergoing an intense and short 
phase of star formation, at a rate significantly higher than the 
average SFR of the galaxy before the starburst event, indepen- 
dently from the galaxy IR luminosity. Adopting Elbaz et al.'s 
(2002 1 definition, we would count 8 LIRGs and 6 'starbursts' 
in our sample of 14 field galaxies. 

The redshift distribution of our field sources (0. 1 ^ z ^ 1 . 1 , 
with a median z - 0.6) is also in agreement with the values 
found in other ISOCAM field surveys (see Flores et al. 119991 
Rigolopolou et al. .2002. Elbaz & Cesarski .2003. Paper I). 

In conclusion, our analysis confirms and reinforces previ- 
ous results, that we can summarize by saying that, on average, 
ISOCAM-selected field galaxies are disk galaxies at z ~ 0.6, 
with ongoing star formation activity at a rate of several dozens 
Mq yr"' . It is possible that some (or even most) of these galax- 
ies have suffered one or more starburst events in the recent past, 
but only a minority seem to be detected during the peak of the 
(short-lived) starburst phase. The contribution of an AGN to 
the IR emission of field galaxies cannot be excluded, but it is 
unlikely to be dominant. 

9. Summary and conclusions 

Using ISOCAM observations of the A2218 cluster (described 
at length in Paper I), optical data from LPS, SOI, and ZOl, and 
HS T imaging, we have investigated the properties of MIR- 
selected sources in the A2218 cluster field. Specifically, we 
have used a battery of GRASIL models (S98) to search for the 
best-fit to the optical-to-MIR SEDs of the non-stellar sources. 
Acceptable and well constrained fits have been obtained for 
the SEDs of 41 sources, of which 27 are cluster members (see 
TableO). 

The properties of the cluster members detected in the MIR 
can be summarized as follows. Most of them are early-type 
galaxies with small or negligible ongoing star-formation, and 
no significant evidence of recent previous episodes of star- 
formation. These galaxies are mostly detected at 6.7 //m only, 
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where we see the Rayleigh-Jeans tail of cold stellar photo- 
spheres. A quarter of the MIR cluster galaxies show significant 
(albeit small) star-formation activity. They are characterized by 
smaller baryonic masses (as indicated by their NIR luminosi- 
ties), on average, and generally have spiral or irregular mor- 
phologies, but two of them are of SO-SOa type. 

The spatial and velocity distributions of the MlR-selected 
cluster galaxies are not significantly different from those of 
the general cluster population. MlR-selected cluster galaxies 
are not, in general, BO galaxies. Several independent analy- 
ses suggest that A2218 is not a dynamically relaxed cluster 
Complementing what is seen in other clusters (see, e.g., C04), 
our results therefore suggest that a cluster active dynamical sta- 
tus is not sufficient to affect the MIR properties of its member 
galaxies. 

Foreground and background galaxies in the A2218 clus- 
ter field are mostly detected at 15 jum. They span the redshift 
range 0.1-1.1, with a median z - 0.6. Their MIR emission 
is most likely dominated by dust-reprocessed stellar radiation, 
although we cannot exclude an AGN contribution in some of 
them. Models of actively star-forming galaxies, such as those 
of normal spirals and starburst galaxies, observed close to the 
time of maximum star-formation activity, but not necessarily 
during the short-lived starburst phase, provide good fits to the 
SEDs of MIR field galaxies. These field galaxies have a wide 
range of morphologies, but most of them are spirals. Features 
suggestive of ongoing or past interaction are found in a minor- 
ity of the galaxies of the sample. About half of the field galaxies 
in our sample are LlRGs, and their high IR luminosities trans- 
late into a median SFR of ~ 20 Mq yr ' . 

This paper is the second in a series devoted to the analysis 
of the data of the ISOCAM gravitational lensing deep survey 
(see Paper 1). In a forthcoming paper of these series (Perez et al. 
12004 ^ we will extend our analysis to the lensing clusters A370 
and A2390, which, together with our analyses of the clusters 
C10024-I-1654 (C04) and A2219 (Coia et al. l2004bi . will pro- 
vide the most extensive data-sample of MlR-selected galaxies 
in medium-distant clusters. This sample will be expanded by 
forthcoming observations with the Spitzer Space Telescope. 
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Appendix A: The SED fitting procedure 

As described in § |3l we fit the SEDs of our ISOCAM sources 
with 30 GRASIL models using a x^ method. Here we detail 
some aspects of the SED fitting procedure. 

The wavelength resolution of the GRASIL models in the IR 
is such that we must take the band- widths of the LW2 and LW3 
ISOCAM filters (which are significantly larger than the band- 
widths of optical and NIR filters) into account, yet, on the other 
hand the detailed shape of the filter wavelength responsivities 



(see, e.g., Blommaert et al. 1200 It is not important, hence we 
consider them to be box-shaped. 

The model SEDs are redshifted to the source z. When the 
source z is not known, we search for the best-fit among 12 val- 
ues of z between z - and z = 3, roughly equally spaced 
in cosmological time, and including the mean cluster redshift, 
z - 0.175. Note that, at variance with Mann et al. J2002> . 
we only consider models with an age less than the age of the 
Universe (in the adopted cosmology) at the given z- 

When the redshift of the optical counterpart is known, there 
are two free parameters in the fit, the model SED and the 
normalisation of the model. In principle, parameters in the 
GRASIL code depend on the galaxy mass (e.g. the duration of 
the wind phase), so that the flux scale cannot be re-normalized 
at will. However, a fine-tuning of the GRASIL parameters is 
only appropriate when the galaxy SED we are trying to fit is 
very well constrained. Otherwise, it is sensible to leave the flux 
scale normalisation as a free parameter in the fit, in so far as 
the observed galaxy masses are not very different from those 
the models were originally designed to fit. The GRASIL mod- 
els we use here had their (mass-related) parameters set to re- 
produce galaxies of L* luminosity or brighter (L. Silva, private 
comm.), hence they are appropriate for a comparison with our 
ISOCAM galaxies. 

When z is not known, it is taken as an additional free pa- 
rameter in the fitting procedure. In order to check how well we 
can determine the photometric redshift, Zpiwt, of our ISOCAM 
galaxies using the GRASIL models and the available optical- 
to-MlR data, we used the ISOCAM sources with known-z, and 
left z as a free parameter in the SED-fitting procedure. The best- 
fit then provides our Zphoi estimate. The photometric and spec- 
troscopic z's are significantly correlated (Kendall's rank corre- 
lation probability is > 99 %), and the root mean square (rms, 
hereafter) difi'erence between the spectroscopic and photomet- 
ric z's is 0. 10. Hence, we can hope to use the SED-fitting proce- 
dure to determine approximate Zpiwi estimates for the ISOCAM 
sources with the best-sampled SEDs. 

In particular, for a given ISOCAM source we consider our 
Zphoi estimate sufficiently well constrained when the rms of the 
z-values obtained in all the acceptable fits (at the > 5 % c.l.) is 
Zrms ^ 0.2. In this case, if Zpiwi ±Zrms is consistent with the clus- 
ter mean redshift, we assign the source to the cluster, otherwise, 
we adopt the best-fit Zphot estimate. 

Source confusion is not a major issue for most of our A22 1 8 
ISOCAM sources. However, when we suspect that another op- 
tical galaxy could contribute to the MIR flux, we should con- 
sider the summed contribution of both optical counterparts. Of 
course, in this case, we must fit two models (and sometimes 
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two redshifts) instead of one. Hence, the number of free param- 
eters sometimes become too large with respect to the number 
of available data, and we cannot constrain the best-fit solution. 
The available data suggest that we try the two-model fit for 
source nos. 33, 45, 61a, 61b, 67. In none of these cases does 
the best-fit improve with respect to the fit with a single optical 
counterpart. 

Appendix B: Notes on individual sources 

We discuss in some detail here the five sources whose SEDs we 
failed to fit with the 30 GRASIL models considered (see §|3- 

ISO^2218.02. This ISOCAM source is detected in both 
the LW2 and LW3 bands. We identify this ISOCAM source 
with the optical counterpart L664, at a redshift z - 0.53, but we 
note that the nearby galaxy L663, at unknown redshift, could 
contribute substantially to the MIR flux, since its r-band flux 
density is 30% of the flux density of L664. Using L664 as the 
only optical counterpart, we do not find acceptable fits to the 
source optical-to-MIR SED. Formally, the best-fit is obtained 
for a Sa5 model, but the fit is poor mostly because the model 
underestimates the B-band flux. It is possible that the MIR flux 
is provided by both L664 and L663. However, we do not have 
enough data to constrain a fit with two independent models for 
the two galaxies, and a free redshift for L663. 

IS0^2218_28. This ISOCAM source is detected in the 
LW2 band only, and its optical counterpart is L391, the clus- 
ter cD, at z = 0.1720. Its spectral type is that of an early-type 
galaxy, according to LPS. The MIR emission of this ISOCAM 
source is very extended, and probably it is not contributed 
solely by the cD (see Fig.[T)- Using L391 as the only optical 
counterpart, this source SED is best fit by an EmlO model, but 
the fit is poor, mainly because it underestimates the K^ band 
flux. Given the very extended emission from this galaxy, and 
the considerable crowding of galaxies in its vicinity, it is per- 
haps not surprising that it is difficult to fit its SED. Consistent 
aperture photometry is needed to ensure that we are sampling 
the same parts of this galaxy at different wavelength. An accu- 
rate modelling of the cD is however beyond the scope of this 
paper 

IS0^2218J5. This is an ISOCAM source detected in 
both the LW2 and LW3 bands. It is identified with L317 (also 
S420), an Scd at z = 0.474. The best-fit to the source SED is 
given by a Mps5 model, but the fit is only marginally accept- 
able (at the 1.6%c.l.). 

IS0^2218^3. This ISOCAM source is detected in the 
LW3 band only. Its most likely optical counterpart is L205 (also 
S368), but L212 might contribute significantly to the MIR flux, 
since the r-band flux ratio of the two galaxies is 0.4. L205 is an 
Sc at z = 0.693, L212 is an unknown-z spiral. Using L205 as 
the only optical counterpart, the best-fit is obtained for a Sc5 
model, but the fit is poor, mainly because it underestimates the 
Ks band flux. Unfortunately, we do not have enough data to 
constrain a fit with two independent models for the two galax- 
ies, and a free redshift for L212. 

IS0^2218.67. This ISOCAM source is detected in both 
the LW2 and LW3 bands. It is identified with L75 (also S4004), 
a faint object at unknown z- Other galaxies could contribute 



to the MIR emission, in particular L71 (also S4005), an Scm, 
but possibly also L67. The r-band flux ratios of L75 to L71 
and to L67 are 1.4 and 0.6, respectively. Adopting L75 as the 
only optical counterpart, we do not find an acceptable fit to the 
source SED, whatever the source redshift. Formally, the best- 
fit is found for an Es2 model at z — 1.5, but it is a poor fit 
that severely underestimates the MIR fluxes. We try adding L71 
as an additional counterpart. Unfortunately, no acceptable fits 
are found in this case either. The formal best-fit is obtained 
for a combination of a z = 1.1 E15 model, and a z = 2.8 
Es2 model. Although the fit now reproduces the LW2 flux, it 
still underestimates the LW3 flux. The role of another possible 
counterpart, L67, cannot be evaluated for lack of photometric 
data. Notably, this source has also been detected in the X-ray 
by Chandra (Sanchez, private comm.). Hence, part of the MIR 
emission could arise from an AGN component embedded in a 
dusty torus, something that GRASIL models do not currently 
account for. 
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